Optimizing global CO emission estimates using a four-dimensional variational data assimilation system and surface network observations

We apply a four-dimensional variational (4D-VAR) data assimilation system to optimize carbon monoxide (CO) emissions for 2003 and 2004 and to reduce the uncertainty of emission estimates from individual sources using the chemistry transport model TM5. The system is designed to assimilate large (satellite) datasets, but in the current study only a limited amount of surface network observations from the National Oceanic and Atmospheric Administration Earth System Research Laboratory (NOAA/ESRL) Global Monitoring Division (GMD) is used to test the 4D-VAR system. By design, the system is capable to adjust the emissions in such a way that the posterior simulation reproduces background CO mixing ratios and large-scale pollution events at background stations. Uncertainty reduction up to 60 % in yearly emissions is observed over well-constrained regions and the inferred emissions compare well with recent studies for 2004. However, with the limited amount of data from the surface network, the system becomes data sparse resulting in a large solution space. Sensitivity studies have shown that model uncertainties (e.g., vertical distribution of biomass burning emissions and the OH field) and the prior inventories used, influence the inferred emission estimates. Also, since the observations only constrain total CO emissions, the 4D-VAR system has difficulties in separating anthropogenic and biogenic sources in particular. The inferred emissions are validated with NOAA aircraft data over North America and the agreement is significantly improved from the prior to posterior simulation. Validation with the Measurements Of Pollution In The Troposphere (MOPITT) instrument version 4 (V4) shows a slight improved agreement over the well-constrained Northern Hemisphere and in the tropics (except for the African continent). However, the model simulation with posterior emissions underestimates MOPITT CO total columns on the remote Southern Hemisphere (SH) by about 10 %. This is caused by a reduction in SH CO sources mainly due to surface stations on the high southern latitudes

Comparing optimized CO emission estimates using MOPITT or NOAA surface network observations

[1] This paper compares two global inversions to estimate carbon monoxide (CO) emissions for 2004. Either surface flask observations from the National Oceanic and Atmospheric Administration Earth System Research Laboratory (NOAA/ESRL) Global Monitoring Division (GMD) or CO total columns from the Measurement of Pollution in the Troposphere (MOPITT) instrument are assimilated in a 4D-Var framework. Inferred emission estimates from the two inversions are consistent over the Northern Hemisphere (NH). For example, both inversions increase anthropogenic CO emissions over Europe (from 46 to 94 Tg CO/yr) and Asia (from 222 to 420 Tg CO/yr). In the Southern Hemisphere (SH), three important findings are reported. First, due to their different vertical sensitivity, the stations-only inversion increases SH biomass burning emissions by 108 Tg CO/yr more than the MOPITT-only inversion. Conversely, the MOPITT-only inversion results in SH natural emissions (mainly CO from oxidation of NMVOCs) that are 185 Tg CO/yr higher compared to the stations-only inversion. Second, MOPITT-only derived biomass burning emissions are reduced with respect to the prior which is in contrast to previous (inverse) modeling studies. Finally, MOPITT derived total emissions are significantly higher for South America and Africa compared to the stations-only inversion. This is likely due to a positive bias in the MOPITT V4 product. This bias is also apparent from validation with surface stations and ground-truth FTIR columns. Our results show that a combined inversion is promising in the NH. However, implementation of a satellite bias correction scheme is essential to combine both observational data sets in the SH